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Coronavirus particles are enveloped and pleomorphic and are thus refractory to crystallization and sym- 
metry-assisted reconstruction. A novel methodology of single-particle image analysis was applied to selected 
virus features to obtain a detailed model of the oligomeric state and spatial relationships among viral 
structural proteins. Two-dimensional images of the S, M, and N structural proteins of severe acute respiratory 
syndrome coronavirus and two other coronaviruses were refined to a resolution of ~4 nm. Proteins near the 
viral membrane were arranged in overlapping lattices surrounding a disordered core. Trimeric glycoprotein 
spikes were in register with four underlying ribonucleoprotein densities. However, the spikes were dispensable 
for ribonucleoprotein lattice formation. The ribonucleoprotein particles displayed coiled shapes when 
released from the viral membrane. Our results contribute to the understanding of the assembly pathway 
used by coronaviruses and other pleomorphic viruses and provide the first detailed view of coronavirus 


ultrastructure. 


Coronaviruses are an important family of human and veter- 
inary pathogens that can cause enteric and respiratory infec- 
tions. Coronavirus infection can lead to respiratory failure, 
gastroenteritis, nephritis, hepatitis, encephalitis, and progres- 
sive demyelinating disease (45, 46). Prior to the emergence of 
severe acute respiratory syndrome coronavirus (SARS-CoV), 
human coronaviruses were primarily viewed as mundane sea- 
sonal pathogens. That perception changed indelibly in 2003 
with the epidemic of severe acute respiratory syndrome that 
emerged from China. 

Coronaviruses are classified with the toroviruses, arterivi- 
ruses, and roniviruses in the order Nidovirales. Nidoviruses 
share some common aspects of RNA synthesis but diverge in 
virion morphology. Arteriviruses are spherical and may contain 
an icosahedral nucleocapsid (81); ronivirus particles are rod- 
like or pleomorphic (12); torovirus particles are toroidal or 
pleomorphic (2, 39); and coronaviruses are typically spherical 
or pleomorphic (16). Coronaviruses derive their name from 
the protruding oligomers of the spike glycoprotein (S), which 
form a coronal fringe around the virion. Nucleoprotein mole- 
cules (N) and the positive-stranded genome comprise the ri- 
bonucleoprotein (RNP) core. The major component of the 
viral membrane is the triple-pass transmembrane matrix gly- 
coprotein (M), which is essential for virus assembly (reviewed 
in reference 45). 

Coronavirus assembly is localized at membranes of the en- 
doplasmic reticulum-Golgi intermediate compartment (42) 
and is mediated by species-specific interactions of the M pro- 
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tein with S, N, and envelope (E) proteins. The stringency of 
structural protein organization at the site of budding is such 
that some transmembrane host proteins resident at the site of 
assembly are excluded from the virion, presumably due to tight 
M-M interactions (18). 

Our general understanding of virion organization has 
been derived primarily from X-ray crystallography and elec- 
tron microscopic analysis of homogeneous, symmetric vi- 
ruses, including viruses with helical or icosahedral capsids, 
and enveloped viruses that display icosahedral symmetry 
(19, 47). Electron micrographs of coronaviruses show par- 
ticles that are neither homogeneous nor symmetric (16, 52, 
64, 67). Thus, a novel image analysis methodology was ap- 
plied to overcome some of the limitations posed by pleo- 
morphic particles. Some techniques applied here were orig- 
inally developed for electron cryomicroscopy (cryo-EM) 
analysis of pleomorphic arenavirus particles (62) and the 
envelope of hepatitis B virus (20). 

Cryo-EM and image analysis were used to examine the su- 
pramolecular structure of three coronaviruses: SARS-CoV, 
feline coronavirus (FCoV), and murine hepatitis virus (MHV). 
Single-particle analysis was applied to groups of images col- 
lected from similar regions of the virion using iterative align- 
ment and averaging (29). Structural information was extracted 
from hundreds to thousands of similar partial virion images. 
Information on vertical organization of viral features was de- 
rived from images boxed at the virion edge, designated “edge 
views,” and complementary information on the horizontal dis- 
tribution and organization of features was derived from “axial 
views,” also known as “en face” views, selected near the virion 
center. In this report, we present a detailed model of the 
supramolecular architecture of coronaviruses based on inte- 
gration of structural data from such axial and edge views. 
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TABLE 1. Number of images selected for processing (x10), 
grouped by virus and orientation 


Virus 
Orientation 
SARS-CoV FCoV MHV TUN-MHV 
Edge 11.0 5.0 3.0 1.5 
Axial 10.0 4.5 2 1:5 


MATERIALS AND METHODS 


Virus growth and purification. Vero-E6, feline AK-D, and murine DBT cells 
were maintained in Dulbecco’s modified Eagle medium (Invitrogen) supple- 
mented with 10% fetal bovine serum (HyClone), 2 mM L-glutamine, and anti- 
biotics. SARS-CoV-Tor2 was propagated on Vero-E6, FCoV-Black was propa- 
gated on AK-D, and MHV-OBLV60 was propagated on DBT cells. Viral 
supernatants were collected 24 h postinoculation, and infectious titers were 
determined by plaque assay as described previously (63). Virions were harvested 
from infected cell culture medium by precipitation with 10% polyethylene glycol 
(PEG-8000). The PEG pellet was resuspended in 0.9% NaCl with 10 mM 
HEPES and applied to a 10 to 30% sucrose density gradient. Final purification 
was achieved by ultracentrifugation at 27,500 rpm for 2 h at 5°C. Virus particles 
resuspended in HEPES-buffered saline were fixed with 10% (SARS-CoV) or 1% 
(FCoV and MHV) formalin overnight at 4°C prior to cryopreservation. 

Cryo-EM. Low-dose cryo-EM was performed at 120 kV as described previ- 
ously (21). The electron dose was minimized to reduce specimen damage. Digital 
cryo-EM images were recorded by two methods. SARS-CoV film negatives were 
digitized using a Zeiss SCAI microdensitometer at a final resolution of 1.84 
A/pixel. All other images were recorded directly via a CCD image sensor at a 
resolution of 2.26 A/pixel using Leginon (75). The signal strength and signal-to- 
noise ratio differed subtly depending on recording methodology. Image contrast 
was inverted so that protein density appeared white. The density histogram for 
each image was normalized to a common median gray value prior to analysis. 

Image analysis. Image processing was performed using the EMAN (53) and 
SPIDER (30) software suites. The EMAN module “boxer” was used to manually 
select regions for analysis (Table 1). Damaged virions were excluded from image 
selection. The EMAN module “ctfit” was used to estimate the contrast transfer 
function parameters for each micrograph. At underfocus conditions, the contrast 
transfer function imposes phase reversals in different spatial frequencies. To 
eliminate this artifact, particle images were subjected to a low-pass Fourier filter 
to truncate frequencies beyond the first phase reversal in the contrast transfer 
function (between 1.6 and 3.5 nm~*). 

The x, y origin and rotational orientation of boxed edge and axial views were 
aligned by 10 rounds of iterative centering and averaging. The EMAN module 
“startnrclasses” grouped images using factor analysis with k-means grouping to 
identify covariant features present in particle images. The extracted covariant 
features were used as a basis for image classification to produce “reference-free” 
class averages. Particles were classified in groups of 50 to 500 edge views or 100 
to 10,000 axial views. Averages for each class were generated by an iterative 
algorithm, where the average from the previous round was used to align particles 
for the next average. For each round of alignment, particles that deviated by 1 
standard deviation (SD) or more from the mean were excluded from the average 
for that particular round. About 60% of aligned input images in a group were 
averaged to produce a final class average. 

Radial density profiles. Although the particles were pleomorphic and not 
circular, 30° sectors within the particles could be identified that had a constant 
radius of curvature. Radial density profiles for individual sectors were calculated 
using X3D (9, 10). The bilayer density was used as a fiducial mark for alignment 
of radial density profiles. 

PCA. Axial images were refined by iterative cycles of alignment. Principal 
component analysis (PCA) was performed to analyze the sources of interimage 
variation. Principal components (eigenvectors) were extracted from the matrix of 
aligned axial images using the SPIDER command CA S. The term eigenvector 
refers to one of a set of ranked factors that describe the prominent covariances 
found in the input matrix [29]. Each eigenvector was converted to visual form, 
known as an eigenimage. Eigenimages were reconstituted from the individual 
eigenvectors using the command CA SRE. The command CA SR was used to 
reconstitute an approximation of the original axial view images using selected 
eigenvectors. The relative contribution of each eigenvector to the reconstituted 
image was weighted in proportion to the prominence of that eigenvector in the 
original image. 
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RESULTS 


Ultrastructure of frozen hydrated coronavirus particles. 
Cryo-EM images of coronavirus particles showed fields of en- 
veloped, spherical, and ellipsoidal virions (Fig. 1 and 2). A 
fringe of cone-shaped spikes protruded from the viral enve- 
lope. The virion core contained ~5-nm punctate densities. 
Particle diameters ranged from 50 to 150 nm, excluding the 
spikes, with mean particle diameters of 82 to 94 nm. 

The precision of other measurements was assessed by mea- 
suring the distance between the inner and outer bilayer head- 
group densities from empty vesicles that were copurified with 
virions. Bilayer headgroup densities were separated by 3.8 nm 
with a standard error of +0.1 nm (0.5 pixel) and an SD of +0.5 
nm, or 2.5 pixels (n = 45). Published studies indicate that the 
electron-dense phospholipid headgroup regions within bilayers 
are separated by 3.6 to 4.0 nm in most common cellular lipid 
membranes (58). Thus, from the coronavirus bilayer measure- 
ment, it was determined that densities in these images could be 
measured with subnanometer accuracy and precision. 

Characterization of viral spikes. Spikes extended ~19 nm 
(SARS-CoV and MHV) to 21 nm (FCoV) from the high- 
density headgroup region of the outer bilayer leaflet. Tunica- 
mycin (TUN) treatment induces misfolding of the spike pro- 
tein (S) and leads to the formation of spike-depleted virions 
(68, 69). TUN-MHV therefore served as a spike-depleted con- 
trol. Spike-depleted TUN-MHV (Fig. 1E) and TUN-SARS- 
CoV particles (not shown) displayed otherwise normal “coro- 
navirus-like” morphology. The mean diameters of TUN-grown 
and native MHV particles did not differ significantly (TUN- 
MHV [means + SD], 93 = 17 nm [n = 87]; MHV, 89 + 13 nm 
[n = 80]; P > 0.05 by ¢ test). Thus, the surface spike is dis- 
pensable for formation of authentic virus-like particles, as pre- 
viously reported (11, 34, 38, 80). 

Characterization of proteins M and N. The appearance of 
closely packed densities in the membrane and the punctate 
electron-dense features in the core of all three coronaviruses 
supported attribution to the remaining, conserved, high-copy 
structural proteins M and N, respectively. M is known to be 
membrane embedded, and the carboxyl tail interacts specifi- 
cally with N (26, 43, 59, 74). The M-N interaction should 
therefore constrain some N molecules in close apposition to 
the envelope. To characterize the M-N interaction, density was 
plotted as a function of radial distance relative to the viral 
membrane. Fewer than 1% of the imaged virions were found 
to be sufficiently circular to permit successful whole-particle 
radial density analysis. However, analysis of the density within 
30° virion sectors yielded more consistent results. Two high- 
density features were consistently detected on the inner bilayer 
leaflet of virion images (Fig. 3A). The 1-nm layer of density 
directly apposed to the inner bilayer headgroup was ascribed to 
M. The second feature, centered 12 nm below the inner head- 
group density, was ascribed to N in the form of the viral RNP. 
The centers of mass for interacting M and N molecules were 
therefore separated by ~10 nm. 

Analysis of structural protein spacing. The Fourier trans- 
formation is a mathematical operation that resolves a signal 
into amplitude and phase components as a function of fre- 
quency (for a review of single-particle EM techniques, includ- 
ing Fourier transformation, see reference 29). Fourier trans- 
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FIG. 1. Cryo-EM of coronaviruses in vitreous ice. SARS-CoV-Tor2 (A and B), FCoV-Black (C), MHV-OBLV60 (D), and tunicamycin-grown 
MHV-OBLV60 (E) are shown in “reversed” contrast, with density in white. Images were recorded at ~2.5 microns below true focus (B, C, 
right-hand portion of panel D, and E), ~4.0 microns under focus (A), and ~20 microns under focus (left-hand portion of panel D). The scale bar 


represents 100 nm. 


formation can be used to analyze the size and spacing of 
component features within a two-dimensional image, such as 
an electron micrograph. Periodicity within an image will pro- 
duce peaks in the power spectrum, a plot of the amplitudes of 
the Fourier transform versus spatial frequency. Power spectra 
were calculated from images of coronavirus particles that were 


recorded at several focus levels; results from two focus levels 
are shown in Fig. 3B. Adjacent regions of background vitrified 
ice were analyzed for each specimen as controls. The power 
spectra of ice-embedded virions displayed two prominent peaks 
of intensity, which were absent in the power spectra from 
images of the background vitrified ice. A sharp peak at 15 
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FIG. 2. Cryo-EM of pleomorphic SARS-CoV virions. Shown is a gallery of images at —2 to —4 microns defocus to illustrate surface and interior 


topography. The scale bar represents 100 nm. 


nm‘ and a broader peak at 5 to 8 nm‘ were observed for 
SARS-CoV (Fig. 3B). Power spectra from images containing 
MHV and FCovV particles displayed similar features (data 
not shown). 

The 15-nm' and 5- to 8-nm ' features were also present in 
power spectra of boxed regions of virion images, including 
axial views (data not shown) and edge views (Fig. 3C). The 
physical separation of S, M, and the RNP in edge views affords 
an opportunity to dissect the origins of prominent peaks in the 
power spectra. SARS-CoV edge views with central spikes were 
aligned so that the membrane was horizontal. Each image was 
then divided into an “inner” half containing the RNP and an 
“outer” half containing spikes (Fig. 3C, inset). Power spectra 
of spike-depleted TUN-MHV particles lacked the 15-nm~+ 
feature present in native SARS-CoV edge views, which con- 
firmed the assignment of this peak to the spikes (data not 
shown). The 5- to 8-nm~* peak was not related to the spikes 
and instead originated primarily from the RNP. Dissection of 
FCoV and MHV edge views yielded the same attribution of 
densities (data not shown). 

The 15-nm ' peak did not match measurements of the spike 
height (~19 to 21 nm) or width (~10 nm). Therefore, inter- 
spike spacing was investigated as a possible source of the 15- 
nm * peak in the power spectra. To accomplish this, aligned 
edge views were duplicated and divided into outer, middle, and 
inner thirds, containing spikes, membrane, and RNP, respec- 
tively (Fig. 3C, inset). One copy of each image was shifted 
horizontally in 1-pixel increments and compared to the second, 
static copy, with the results expressed as a correlation coeffi- 
cient. Correlation coefficients are expressed on a scale of —1 to 
1, with positive values indicating positive correlation. Correla- 
tion peaks at 13 to 15 nm in the outer and middle regions 
reflected the most commonly observed interspike spacing (Fig. 
3D). A small increase in correlation between 5 and 8 nm in the 


RNP region was insufficient to characterize the inter-RNP 
spacing. 

Lattices of spike and RNP densities. Analysis of the edge 
views revealed that the 15-nm~' spike and 5- to 8-nm~' RNP 
features in the power spectra were not interdependent, so they 
could be analyzed separately (Fig. 3C and D). Analysis of axial 
views of spike-depleted and native coronavirus particles pro- 
vided an independent approach for examining the structure of 
the envelope without any assumptions regarding the symmetry 
or packing of the glycoprotein surface spikes. By extracting 
areas in the center of the particle, the radius of curvature is 
sufficiently small that, at least at low resolution, we can treat 
these areas as if they were planar. Axial images represent the 
superposition of the near and far sides of the virus as well as 
any interior density. The process of reference-free alignment 
was presumed to favor the less distorted or more prominent 
side of the virion, although we were not able to test this 
assumption directly. Fourier space filtration was therefore 
applied to resolve the spike and RNP contributions in axial 
views, where the two features overlap (Fig. 4). Axial images 
were filtered to remove image components either larger or 
smaller than 9 nm. Filtered images were then iteratively 
aligned and averaged in real space. Images of background 
ice were processed similarly to control for artifacts that 
might have been introduced during the filtration process. 
Features larger than 9 nm were refined as round, 10-nm- 
diameter spikes. The spike spacing approximated an oblique 
planar crystallographic lattice with a unit cell with edges of 
~14 and 15 nm and an angle of ~100° (Fig. 4C, row c). 
Analysis of axial images of spike-depleted TUN-MHV 
showed no lattice, confirming the assignment to the glyco- 
protein spikes (Fig. 4C, column 4, row c). 

Features smaller than 9 nm in SARS-CoV, FCoV, MHV, 
and TUN-MHV axial images were refined as arrays of oval 
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FIG. 3. Disposition of densities in the membrane-proximal region. (A) Rotationally averaged radial density profiles were generated for 30° 
wedges taken from intact coronavirus particles. Wedges from SARS-CoV (n = 80), FCoV (n = 41), MHV (n = 53), and TUN-MHV (n = 82) 
particles were aligned on the minimum-density node between the headgroup densities of the lipid bilayer. The schematic above interprets densities 
in the spike protein (S), membrane-proximal and matrix protein (M), and ribonucleoprotein (RNP) regions. Entire SARS-CoV virions and 
adjacent regions of background vitrified ice, recorded at two levels of focus, are presented as reciprocal space power spectra (B). Prominent 
features are noted at ~15 nm“! (heavy red bracket) and ~5 to 8 nm‘ (dotted bracket). An averaged power spectrum for ~1,000 aligned 
SARS-CoV edge views (C) or the extraviral (outer) or interior (inner) portion of the image is also shown. The scale bar represents 10 nm. 
Prominent features are marked as described for panel B. Edge views (C, entire) were divided into upper (D, outer), membrane-proximal (D, 
middle), or lower (D, inner) thirds. The correlation coefficient was calculated for a static copy of each image and another copy of the same image 
that was shifted in 1-pixel increments. Peaks of correlation are marked as described for panel B. 


RNP densities. RNP densities of 5 by 6 nm were arranged in an 
~100° oblique lattice, with unit cell edges of 6 and 7.5 nm (Fig. 
4C, row a). Fourier transformation of the spike and RNP 
lattices yielded first-order reflections (insets), consistent with 
locally ordered planar packing. Similar reflections also origi- 
nated from the centers of some individual particles (Fig. 4D). 
The 6.0- to 7.5-nm inter-RNP reflections imply medium-range 
“paracrystalline” packing, with consistent spacing but variable 
alignment of fine features. The RNP densities were consistent 
with the dimensions and appearance of coronavirus RNP re- 
leased from disrupted virions (15, 55). 

We then analyzed the superposition of lattices of the spike 
and RNP features. Reflections from Fourier-transformed class 
averages of axial views were isolated and used to reconstruct 
the spike and RNP lattices. Reconstructed SARS-CoV, FCoV, 


and MHV spikes and RNP lattices were found to be nearly 
identical. Each spike was positioned centrally over four RNP 
densities (Fig. 5). 

The resolution limit of the axial class average images was 
next examined. Class averages were refined using only the 
even- or odd-numbered images in a class and then compared 
by Fourier ring correlation (FRC) analysis (71). The point 
where the FRC value equaled 0.5 was taken as the nominal 
resolution limit. Arrays of RNP densities in axial class averages 
were resolved to 5.3 to 5.8 nm (Fig. 4C and D and 5A). 

Oligomerization of S. Principal component analysis (PCA) is 
a statistical technique that can be used to examine interimage 
variance. Here, PCA was applied to improve the resolution of 
spikes in the axial images. Coronavirus spikes have been re- 
ported to be homotrimeric based on structures of the fusion- 
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FIG. 4. Analysis of viral features in axial views. (A) Axial spike 
images were selected from the central region of each virion (left 
image) and masked around one spike (center and right images). 
(B) An averaged image from unaligned, unmasked, manually centered 
SARS-CoV axial views is also shown. (C) Axial images of SARS-CoV 
(column 1), FCoV (column 2), MHV (column 3), and TUN-MHV 
(column 4) were aligned and averaged iteratively until a stable aver- 
aged image emerged (row b). Axial images were filtered in Fourier 
space to remove image data greater than (row a) or less than (row c) 
9 nm. Filtered axial images were averaged, the averaged image was 
refined by 10 rounds of iterative alignment and averaging, and then 
unfiltered images were aligned to the averaged filtered image for a 
further two cycles to produce the images shown (rows a and c). Insets 
show Fourier transforms of the corresponding averaged images. 
(D) One virion that displayed reflections corresponding to the spike 
and RNP lattices is shown, together with a Fourier transform of the 
core region (inset). The scale bar represents 10 nm. 


activated core (24, 76, 79, 85, 86). However, some reports have 
described homodimeric forms of the spike (7, 48). 

Some reference-free class averages of axial views appeared 
to show triangular spikes (Fig. 5E, left). To investigate further, 
a SARS-CoV reference-free class average was refined using 
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~1,000 axial images. Principal components (eigenvectors) 
were extracted from the matrix of aligned images. 

Prominent eigenimages showed ~10-nm-diameter spikes, 
situated ~15 nm apart (one example is presented in Fig. 5E, 
center). Individual axial-view images were reconstructed from 
the weighted contributions of four prominent eigenvectors 
showing spikes. Three-lobed spike densities were recon- 
structed from axial SARS-CoV images (Fig. SE, right). In- 
tralobe spacing within the trimer was ~6 nm. Reconstructed 
FCoV and MHV spikes also showed similar three-lobed den- 
sities (data not shown). 

Images of ice and axial TUN-MHV images were analyzed by 
PCA as a control. Eigenimages from background ice resem- 
bled featureless noise. Axial TUN-MHV eigenimages and re- 
constructed images showed only the RNP lattice (Fig. 5F). The 
observed three-lobed densities therefore originated only from 
images showing spikes. 

The approximate volume occupied by a spike was estimated 
as another way to assess the oligomeric state of the S protein 
in the spikes. The average partial specific volume of folded 
protein, measured at ~0.73 cm?/g (36), was used to calculate 
the expected volume of the spike ectodomain. One copy of a 
SARS-CoV (130 kDa), FCoV (151 kDa), or MHV (137 kDa) 
S protein ectodomain was therefore predicted to occupy 1.8 X 
10? to 2.1 x 10° nm?, for a hypothetical volume of 5.4 x 10? to 
6.3 X 107 nm? per trimer. Axial and edge projections of the 
spike were measured from class averages and images that had 
been filtered to remove high-resolution noise. On the basis of 
the dimensions of the spikes, the volume of an ellipsoidal spike 
“head” on a cylindrical “stalk” was estimated at 4.9 x 10° to 
5.9 X 10° nm°. Spikes that were modeled instead as cone 
segments gave estimated volumes of 5.2 X 10% to 5.8 x 10? 
nm. The volume of each spike ectodomain is therefore con- 
sistent with a homotrimer. 

Estimation of S:N stoichiometry. The volume of the RNP 
densities was next estimated. RNP densities modeled as ellip- 
soids gave estimated volumes of 53 to 73 nm*. The predicted 
protein volume for the N monomer of SARS-CoV, FCoV, or 
MHV is 51 to 60 nm? (36). Each RNP density therefore con- 
tains one N protein molecule. The native coronavirus spike was 
thus found to be trimeric and was situated over four RNP 
densities, each likely consisting of one nucleoprotein density 
plus RNA. The ratio of S to N at the virion surface was 
therefore found to be three copies of S to four RNP densities 
(S3:4N). The observed structural unit of 1 xX S,:4N is also 
consistent with the dimensions of the spike and RNP unit cells. 

Vertical organization of the S-M-RNP complex. Edge view 
images were refined to examine the stoichiometry of S:M and 
M:N. M is present in higher copy numbers than N or S in the 
virion (6, 25, 35), and it likely plays a pivotal role in assembly. 
Well-defined edge-view class averages were obtained by mask- 
ing out parts of the images that showed the highest variability. 
Empirical testing revealed that a keystone-shaped mask (Fig. 
6A and B) incorporating one spike and two adjacent RNP 
densities of the membrane-proximal shell (Fig. 6C) gave the 
optimal results. The area refined most clearly corresponded to 
a side projection of the unit cell that was identified in axial 
views. 

Refined FCoV edge views (Fig. 6D) showed spikes with 
elongated stalk regions that may correspond to the extended 
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FIG. 5. Analysis of SARS-CoV paracrystalline lattices. The SARS-CoV RNP lattice (Fig. 4C, column 1, row a) was used as the starting image 
for 10 cycles of alignment and averaging (A, left image). Reflections were selected from the Fourier transform of this image (A, right) and back 
transformed to reveal the overlapping RNP (B, left) and spike (B, right) lattices. (C) Lattices similarly obtained from reflections in Fourier 
transformed FCoV (left) and MHV images (right) show similar RNP organization. (D) The two lattices are shown superimposed and averaged 
(left) and in schematic form (right). (E) A reference-free class average of axial SARS-CoV images recorded slightly farther from focus shows 
triangular spikes (left), also seen in the example principal component eigenimage (center) and in reconstructed images (right). (F) An example 
eigenimage from TUN-MHV, showing only RNP densities, is presented for comparison. The scale bar represents 10 nm. 


amphipathic heptad-repeat regions predicted to reside within 
the stalk region (4). Refined edge views of TUN-MHV parti- 
cles showed multiple RNP densities with no external spikes 
(Fig. 6F). Spikes were resolved most clearly in images with only 
one layer of resolved RNP densities (Fig. 6C to E). Analysis of 
interimage variability confirmed that the spikes are in register 
with only the RNP densities in the outermost shell and that 
regions outside the aligned densities are highly variable (data 
not shown). 

The resolution of edge-view class averages was measured by 
FRC analysis as described for the axial views. The nominal 
resolution of the SARS-CoV and FCoV edge class averages 
was 4.2 to 4.7 nm. TUN-MHV edge-view class averages were 
resolved to 3.8 to 4.3 nm. MHV edge-view class averages were 
more poorly resolved (>6 nm), likely due to glycoprotein shed- 
ding during purification. 

Characterization of M. Intramembrane densities were 
poorly resolved in class averages. However, intramembrane 
densities spaced 3 to 5 nm apart were observed in some raw 
virion images (Fig. 6G and H). Similar intramembrane densi- 
ties were observed in spike-depleted TUN-MHV. M is known 
to be a triple-pass integral membrane protein (18). Therefore, 
the intramembrane densities were related to transmembrane 
regions of M rather than S. Densities spaced 5 to 8 nm apart 
connected the RNP to the membrane region (Fig. 6G and H, 
white arrowheads), consistent with the spacing of paracrystal- 
line RNP. The 5- to 8-nm connecting densities could not be 
ascribed confidently to M or N. 

Spherical and coiled RNP configurations. Edge views re- 
vealed that most of the viral RNP was located within ~25 nm 
of the inner face of the membrane (Fig. 6A). The diameter of 
coronavirus RNP complexes released from disrupted particles 
has been reported to be 14 to 16 nm (15, 55). Fifteen-nano- 
meter-wide strands of electron-dense material can be seen 
emerging from a spontaneously disrupted SARS-CoV particle 
(Fig. 61). The appearance of the 15-nm strands is consistent 
with published descriptions of helical-form coronavirus RNP 


(15, 55). The RNP is maintained in a spherically packaged 
form at the inner face of the membrane (Fig. 61, large and 
small black arrowheads). Flat-edged cores, three- and fivefold 
symmetry axes, and other hallmarks of icosahedral symmetry 
were not detected in the images of frozen-hydrated coronavi- 
rus particles. These results did not rule out the possibility that 
local paracrystalline ordering could give rise to the cores with 
some planar surfaces reported by others (67). 


DISCUSSION 


In this study, cryo-EM and image analysis were used to 
examine the structure of three coronaviruses: SARS-CoV, 
FCoV, and MHV. Our results demonstrate a higher level of 
supramolecular organization in the viral envelope region than 
was evident from previous EM studies, and they therefore 
suggest some modifications to previous interpretations of coro- 
navirus structure (23, 54, 73), as shown in Fig. 7. Glycoprotein 
spikes were found to be aligned with the membrane-proximal 
layer of RNP densities, implying that protein location within 
the envelope is constrained by consistent S-M, M-M, and M-N 
interactions. The clustering of viral features near the mem- 
brane was consistent with published biochemical data showing 
S-M (17, 34), M-M (18), and M-N (26, 43, 59, 74) interactions. 
Although some particles displayed membrane densities ascrib- 
able to M (Fig. 6G and H), the relative visualization of M in 
cryo-EM images of whole virions may have been limited by the 
signal-to-noise ratio (Fig. 3B). Visualization of M would also 
be obscured by the large signal arising from the phospholipid 
headgroups of the viral membrane in edge views. However, M 
interacts with both S and N. Results presented here confirm 
the observation of Risco and collaborators (67) that the coro- 
navirus nucleocapsid is separated from the envelope by a gap, 
which we have revealed to contain thread-like densities that 
connect the M protein density on the inner face of the viral 
membrane to a two-dimensionally ordered ribonucleoprotein 
layer (Fig. 6G and H). Organization of S and N in related 
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FIG. 6. Analysis of the structural proteins as seen in edge views. (A) Boxed images were centered on the viral membrane below one spike, as 
shown in the averaged image of aligned SARS-CoV edge views. (B) A keystone-shaped mask (shown in panel A) was applied to individual edge 
views before reference-free classification and averaging. (C) Edge-view class averages show the ultrastructure of the membrane-associated 
structural protein complex. Class averages of FCoV (D), MHV (EB), and spike-depleted TUN-MHV (F) were obtained by the same process. 
Intramembrane densities ascribed to viral M protein are indicated with black arrowheads in SARS-CoV (G) and TUN-MHV (H). Connecting 
densities between the RNP and membrane regions are indicated with small white arrowheads. (I) A spontaneously disrupted SARS-CoV particle 
shows RNP strands escaping from a spherical core, recorded at —4 wm (left) and —2.4 wm (right) under focus. The narrowest parts of the RNP 
are indicated with white arrowheads, and the membrane section remaining with the core is indicated with large black arrowheads (left). The bilayer 
densities are more visible in the near-focus image, indicated with small black arrowheads (right). The scale bars represent 10 nm. 


lattices implies that M is also organized in a two-dimensional 
lattice, and it may provide a scaffold for viral assembly. M 
contains three transmembrane helices (18). The observed ~2- 
nm-wide intramembrane M densities could accommodate a 
bundle containing the three ~1-nm transmembrane helices 
originating from one M protein molecule. However, further 
investigation will be needed to determine the structural basis 
for the inferred M protein scaffold. On the basis of our current 
analysis, we estimate the stoichiometry of the unit cell at the 
virion surface to be approximately 1S,:16M:4N to 183:25M:4N 
proteins, with the remainder of the N protein distributed 
throughout the virion core. 

The presence of overlapping scaffolds of viral proteins in the 


envelope is likely related to the mechanism of viral budding, 
but this does not specifically explain the formation of closed 
spherical virions. Particle closure would require some mecha- 
nism of disrupting the two-dimensionally ordered S:M:N scaf- 
fold. Our analysis did not reveal any consistent structural per- 
turbations that could explain the mechanism of particle 
closure. The single-particle analysis methods applied in this 
study rely on averaging to determine the most common ar- 
rangement of proteins. As such, these methods are well suited 
to analysis of the “guiding principles” of coronavirus particle 
formation but are poorly suited to the analysis of structural 
microheterogeneity. A possible mechanism by which minor 
membrane constituents, such as E and the group-specific inte- 
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FIG. 7. Interpretation of coronavirus virion structure. Conserved 
structural proteins are drawn as they appear in edge views (left) and 
axial views (right). Trimeric spikes (S) are shaded in red, membrane 
proteins (M) are in solid blue, and nucleoproteins (N) are shaded in 
violet. The dimensions of lattices of S trimers (a = 14.0 nm, b = 15.0 
nm, y = 100°) and RNP molecules (c = 6.0 nm, d = 7.5 nm, ¢ = 100°) 
were determined from the reflections shown in Fig. SA and were 
consistent with real-space measurements of the same parameters. 


gral membrane structural proteins, may promote membrane 
curvature is implied by data indicating that mutations in E 
result in a higher proportion of aberrantly budded virions (27). 
Expression of SARS-CoV E was recently shown to induce 
membrane permeabilization (51), and other coronavirus E 
proteins induce the formation of convoluted membranous 
structures similar to those seen in coronavirus-infected cells 
(14, 65). We therefore hypothesize that E might disrupt sheets 
of paracrystalline S, M, and N by intercalating among M pro- 
teins during assembly. MHV lacking E is viable, although 
highly attenuated (44), suggesting that the role of E in viro- 
genesis is either dispensable or that the role of E during bud- 
ding can be compensated by other viral and cellular proteins. 
SARS-CoV particles likely contain at least the three-pass in- 
tegral membrane protein Sars-3a (78, 83) and the type I trans- 
membrane protein Sars-7a (78), in addition to the E protein 
homolog encoded by subgenomic RNA 4. Further investiga- 
tion will be required to determine how the E protein and the 
other coronavirus group-specific proteins contribute to particle 
formation. 

A coronavirus particle of average diameter with a full com- 
plement of spikes could contain ~50 to 100 spike trimers and 
~200 to 400 copies of N in the membrane-proximal lattice 
region, depending on whether the spike or RNP spacing is 
used to calculate the surface area of a spike unit cell. The ratio 
of S:N in surface lattices may decrease if some particles bud 
with less than the full complement of spikes, as in the case of 
TUN treatment. The observation of an oblique lattice of tri- 
meric spikes at the virion surface was unprecedented. Simple 
hexagonal close packing of homotrimeric S should define a 
120° unit cell with a lattice spacing reflecting the 10-nm spike 
diameter. Retroviral particles, for example, appear to contain 
hexagonally arranged proteins (32, 61, 88). Placement of coro- 
navirus spikes is therefore likely controlled by interaction with 
other ordered components of the virion. S is captured at the 
site of budding by the presence of M; thus, interaction with M 
is a likely source of S organization. Nucleoprotein molecules in 
the paracrystalline RNP shell abut, and may be partially orga- 
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nized through, interactions at points of contact in the RNP 
lattice (Fig. 5B and C). The distribution of density in the viral 
core was consistent with a membrane-proximal RNP lattice 
formed by local approaches of the coiled ribonucleoprotein. 
Dimeric N-N interactions (8, 37, 49, 77) are another possible 
source of virion organization. For instance, it was recently 
observed that the N protein dimerization domains of arterivi- 
ruses and coronaviruses share a common fold despite low ho- 
mology (8). Further experimentation will be required to deter- 
mine how N-N interactions contribute to arterivirus core 
morphology. 

Emerging evidence suggests that membrane-proximal pro- 
tein lattices may be characteristic of enveloped non-icosahe- 
dral viruses and may represent an alternative to icosahedral 
and helical symmetry in virion architecture. For example, 
cryo-EM and cryoelectron tomography have demonstrated 
that glycoprotein spikes in the poxvirus core envelope (13, 22) 
and the envelope of hepatitis B virus (20) are hexagonally 
packed. Retroviral Env may be organized similarly: hexagonal 
arrays were described for human foamy virus (82), and loose 
hexagonal packing was noted for human immunodeficiency 
virus type 1 and simian immunodeficiency virus (89). However, 
a subsequent cryoelectron tomography study did not detect 
ordered packing of Moloney murine leukemia virus Env (28). 
In retrovirus particles, Env distribution may be constrained by 
interaction with underlying hexagonal lattices of Gag that have 
been observed for several retroviruses (5, 32, 33, 41, 50, 60, 90). 
Cryo-EM and image analysis of arenavirus particles revealed 
ribonucleoprotein densities that appeared to be organized in a 
two-dimensional, membrane-proximal lattice (62). Radial den- 
sity profiles of influenza A virus particles showed that the 
protein density is highest in the membrane-proximal region, 
perhaps indicative of ordered packing (31). Freeze-fractured 
Sendai virus images showed that the viral matrix protein, which 
promotes viral assembly, is organized in two-dimensional or- 
thogonal arrays at the fracture face of the budding membrane 
(1). Identical lattices of M protein have also been observed at 
the fracture face of the newly formed Sendai virus particles (1). 

Previous reports have noted that surface protein lattices on 
enveloped virions become disordered with increasing time and 
temperature. For instance, lattices of Sendai M protein disap- 
peared upon prolonged storage at 37°C (1), while deterioration 
of vaccinia virus core envelope lattices occurred more slowly 
during cold storage (22). The kinetics of coronavirus lattice 
disruption and effects on infectivity are unknown and will re- 
quire further study. Nevertheless, it appears that viral assembly 
is guided by transient lattice formation, followed by thermal 
disruption of the lattice, rather than a model of nonspecific 
protein packing induced by low temperature. Rapid prepara- 
tion, low-temperature storage, and formaldehyde fixation may 
have contributed to the preservation of paracrystalline lattices 
observed in coronavirus particles. 

The lattices observed on viral particles are therefore best 
interpreted as remnants of the viral assembly process. Infec- 
tious coronavirus particles are assembled upon the arrival of 
the RNP to the pre-Golgi membranes. Although S proteins of 
some coronaviruses are cleaved after budding, cleavage is not 
required per se for infectivity (3), and S proteins of several 
coronaviruses, including SARS-CoV, are predominantly un- 
cleaved (72, 87). Since coronaviruses assemble at pre-Golgi 
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membranes, glycans on the exposed spike ectodomains un- 
dergo additional trimming and modification. The glycosylation 
inhibitor TUN and inhibitors of glucosidase I and II reduce 
coronavirus infectivity, although inhibition of Golgi-resident 
mannosidases does not (66). These results are consistent with 
observations from other pleomorphic viruses, in which addi- 
tion of high-mannose chains and initial glucosidase activity are 
more important than later mannosidase activity in determining 
infectivity (57, 84). The cores of murine and porcine corona- 
virus particles undergo an internal condensation in the trans- 
Golgi that does not coincide with known physical alterations to 
viral proteins, but they may augment infectivity (40, 70). Anal- 
ysis of the cryo-EM images did not reveal any internal features 
within the ~25-nm-thick RNP zone proximal to the envelope. 
This suggests that inner core densities of mature coronaviruses 
are not consistently ordered with respect to the membrane. 
Coronavirus particles appear to bud from discrete patches of 
RNP subjacent to the intermediate compartment membrane 
(56). Results presented in this study suggest that structural 
protein lattice formation is integral to coronavirus budding. As a 
prerequisite to defining the assembly mechanism of pleomorphic 
coronavirus particles, it will be essential to come to an under- 
standing of the assembly pathway that guides lattice formation. 
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